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Abstract 

i . We examine a simple, fuel-air, model of combustion in spark 

| ignition (si) engine with indirect injection. In our two fluid model, 

variations of fuel mass burned in cycle sequences appear due to 
. stochastic fluctuations of a fuel feed amount. We have shown 

| that a small amplitude of these fluctuations affects considerably 

the stability of a combustion process strongly depending on the 
• quality of air-fulel mixture. The largest influence was found in the 

"^5 [ limit of a lean combustion. The possible effect of nonlinearities in 

the combustion process were also discussed. 

c '■ 
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1 Introduction 

Cyclic combustion variability, found in 19th century by Clerk 
(1886) in all spark ignition (si) engines, has attracted great in- 
terest of researchers during last years (Heywood 1988, Hu 1996, 
Daw et al. 1996, 1998, 2000, Wendeker et al. 2003, 2004). Its 
elimination would give 10% increase in the power output of the 
engine (Heywood 1988). The key source of their existence may 
be associated with either stochastic disturbances (Roberts et al. 
1997, Wendeker et al. 1999) or nonlinear dynamics (Daw et al 
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1996, 1998) of the combustion process. Daw et al. (1996, 1998) 
and more recently Wendeker et al. (2003, 2004) have done the 
nonlinear analysis of the experimental data of such a process. 
Changing an advance spark angle they observed the considerable 
increase of the noise level (Wendeker et al. 2003) claiming that 
it is due to chaotic dynamics of the process. On the other hand 
the main sources of cyclic variability were classified by Heywood 
(1988) as the aerodynamics in the cylinder during combustion, 
the amount of fuel, air and recycled exhaust gases supplied to the 
cylinder and a mixture composition near the spark plug. In this 
paper we will model the variation of fuel ignition amount as the 
most common source of instability in indirect injection. 

The present paper is organized as follows. After the intro- 
duction in the present section we define the model by a set of 
difference equations in the next section (Sec. 2). This model, 
in deterministic and stochastic forms, will be applied in Sec. 3, 
where we analyze the oscillations of burned mass. Finally we de- 
rive conclusions and remarks in Sec. 4. 



2 Two fluid model of fuel-air mix- 
ture combustion 

Starting from fuel-air mixture we define the time evolution of the 
corresponding amounts. Namely, we will follow the time histories 
of the masses of fuel rrif, and air m a . 

Firstly, we assume the initial value of m a (i), nif(i) automati- 
cally their ratio r(i): 

m a (i) 
m/(t) 

for i = 0. 

Secondly, depending on parameter r with reference to a sto- 
ichiometric constant s we have two possible cases: fuel and air 
deficit, respectively. For a deterministic model, the first case lead 
to 

r(i) > 1/s (2) 
we calculate next masses using following difference equations: 

m/(i + l) = (1 - a) (mf(i) — ^m a (i)j + Srrif 

m a (i + l) = 5m f, (3) 
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Table 1: Constants and variables of the model. 



stoichiometric coefficient 


s = 14.63 


exhaust ratio 


a = 0.92 


air mass in a cylinder 


m a 


fuel mass in a cylinder 


m f 


fresh air amount 


5m a 


fresh fuel amount 


Sm f 


air/fuel ratio 


r = m a /m,f 


burned fuel mass 


Arrif 


combusted air mass 


Am a 


air /fuel equivalence ratio 


A 


random number generator 


N(0,l,i) 


mean value 




of fresh fuel amount 


5m fo 


standard deviation 


of fresh fuel amount 


cr m f 


standard deviation 




of the equivalence ratio 


OA 
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Figure 1: The combustion curve Am(A) for the constant fresh air feed 
5m a = 200 mg. 



where a is the exhaust ratio of the engine, 5m f and 5m a denotes 
fresh fuel and air amounts added in each combustion cycle i. In 
the opposite (to Eq. 2.2) case 

r(i) < l/s (4) 

we use the different formula 

rrif(i + 1) = 5m f 

m a (i + l) = (1 - a) (m a (i) - smj(i)) + 5m a (5) 

Note that variables m a and mj are the minimal set of our 
interest. From the above equations one can easily calculate other 
interesting quantities as the combusted masses of fuel Amj(i) and 
air Am a (i) and air- fuel equivalence ratio before each combustion 
event i: 

_ ^ m a (i) + (1 - a)Affl„(t - 1) 
~ S mf(i) + (1 - a)Am f (i - 1) 

Basing on experimental results we use the additional necessary 
condition (Kowalewicz 1984) of combustion process 

0.6 < A < 1.3. (7) 
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Figure 2: The dependence of A on sequential cycles i for deterministic 
(a) and stochastic (b) processes. The fresh air amount was assumed to 
be 8m a = 200 mg while the fresh fuel amount varies: 8m f = 13.50 mg , 
14.63 mg, 21.00 mg starting from the top curve, respectively. 



For better clarity our notations of system parameters: constants 
and variables are summarized in Tab. 1. 

Basing on the relations (Eqs. 2.1-2.7) we plotted the com- 
bustion curve for the assumed constant fresh air feed Sm a = 200 
mg. 

Finally, in the case of stochastic injection, instead of constant 
6rrif(i) (Eqs. 2.3 and 2.5) (for each cycle i) we introduce its mean 
value 5m f Q = const., while 8m t in the following way: 

5m f (i) = 5m fo + a mf N(0,l,i), (8) 

where N(0, 1, i) represents random number generator giving a se- 
quence i of numbers with a unit-standard deviation of normal 
(Gaussian) distribution and the nodal mean. The scaling factor 
a m f corresponds to the mean standard deviation of the fuel in- 
jection amount. The cyclic variation of 5mf(i) can be associated 
with such phenomena as fuel vaporization and fuel- injector vari- 
ations. 
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3 Oscillations of burned fuel mass 



Here we describe the results of simulations. Using Eqs. 2.1-2.8 
we have performed recursive calculations for deterministic and 
stochastic conditions and obtained time histories of various system 
parameters: rrif, m a , Amj, Am a and A. The results for A are 
shown in Fig. 2. The upper panel (Fig. 2a), corresponding to 
deterministic combustion for three different values of fuel injection 
parameter 5m f, shows A as straight lines versus cycle i, while 
the lower (Fig. 2b) one reflects the variations of A in stochastic 
conditions. The order of curves appearing in the Fig. 2b is the 
same as in Fig. 2a stating from the smallest value of considered 
fuel injection amounts from the top. In stochastic simulations we 
used input of random fuel injection with standard deviation equal 
to 10% of its mean value a m f = 0.1 5m /<>. The obtained results 
clearly indicate that the fluctuations of A are growing with larger 
A. This can be also found by analytical evaluation of Eq. 2.6. It 
is not difficult to check that 

<?\ ~ A 2 <7 m/ . (9) 

The results for burned fuel mass Am/ are presented in Fig. 3. 
Starting from deterministic conditions (5mf = 5mf = const.) we 
obtain the constant fraction of the burned fuel mass Am/ repre- 
sented by the three straight lines in Fig. 3a lying very close to 
each other. In Fig 3 b-c we show the same, Am/, for the con- 
sidered case of assumed fuel injection {5m / = 13.50 mg - Fig. 
3b, 5m f a = 14.63 mg - Fig. 3c, <5m/ G = 21.00 mg - Fig. 3d) 
and stochastic conditions. Due to different magnitudes parame- 
ter A fluctuations, and dependence of combustion curve Fig. 1 it 
is not surprise that the fluctuations of Am/ have different char- 
acter in all these cases. For lean combustion, which is a stable 
process in deterministic case, the fuel injection fluctuations intro- 
duce considerable instabilities to the combustion process leading 
to the suppression of combustion because in some cycles (Fig. 3b) 
where A is larger that 1.3. Then Equation 2.7 is not satisfied. In 
the next case (Fig. 3c) the effect of stochasticity is much smaller. 
Here we have optimal air-fuel mixture. First of all one should 
note that fluctuations of A are smaller than in previous case (Fig. 
2b). Moreover A oscillate around the region (A ~ 1) in combus- 
tion curve (Fig. 1) which does not have big changes comparing 
to previous case. Finally, Fig. 3d shows the sequence of Am/ for 
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Figure 3: The dependence of burned fuel mass on sequential cycles i for 
deterministic (a) and stochastic (b-d) processes. 5m a = 200 mg while 
8m f takes different values: 13.50 mg , 14.63 mg, 21.00 mg denoted in 
particular figures a-c. 
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the large 5m f (rich fuel-air mixture). The fluctuations of A are 
the smallest of all three ones but A ~ 0.7 causes suppressions of 
combustions in some cycles similarly to the case shown in Fig. 3b. 

4 Conclusions 

In this paper we examined the origin of combusted mass fluctua- 
tions. In case of stochastic conditions we have shown that depend- 
ing on the quality of fuel-air mixture the final effect is different. 
The worse situation is for lean combustion. The consequences of 
it can be observed for idle speed regime of engine work. Unstable 
engine work, interrupted by the cycles without combustion lead 
to a large increase of fuel use. 

Although the presented two component model is very simple 
it can reflect the underlying nature of engine working conditions. 
In spite of fact that the model is characterized by the nonlin- 
ear transform (Eqs. 2.3, 2.5 and 2.7) similar to logistic one, we 
have not found any chaotic region. Possibly that such solutions 
can be found for non realistic model parameters like A and a. 
The other strong limitation was concerned with the sharp edges 
of combustion curve (Amj versus A Fig.l) modelled by a Heav- 
iside step 0(x) function Eq. 2.7. We used such an approxima- 
tion as a simplest one but modeling with the exponential growth 
exp(— 1/x) is more realistic and possible. Similar assumptions of 
the exponential dependence led to chaotic behaviour in papers 
(Daw et al. 1996, 1998, Wendeker 2003). From a physical point 
of view mixture gasoline-air is not uniform before ignition and 
that can cause nonuniform combustion smearing the edges of the 
combustion curve Fig. 1. Calculations considering this effect are 
in progress and the results will be reported in a separate future 
publication. 
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